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The ^^Cl(p,7)'^''Ar reaction rate is of potential importance in the rp-process powering type I X-ray 
bursts. Recently Clement et al. 1] presented new data on excitation energies for low lying proton 
unbound states in ^^Ar obtained with a new method developed at the National Superconducting 
Cyclotron Laboratory. We use their data, together with a direct capture model and a USD shell 
model calculation to derive a new reaction rate for use in astrophysical model calculations. In 
particular, we take into account capture on the first excited state in ^^Cl, and also present a realistic 
estimate of the remaining uncertainties. We find that the ^^Cl(p,7)''^Ar reaction rate is dominated 
entirely by capture on the first excited state in ^'^Cl over the whole temperature range relevant in 
X-ray bursts. In the temperature range from 0.2 to I GK the rate is up to a factor of 70 larger than 
the previously recommended rate based on shell model calculations only. The uncertainty is now 
reduced from up to a factor of 1000 to a factor of 3 at 0.3-0.7 GK and a factor of 6 at 1.5 GK. 



I. INTRODUCTION 



Proton capture rates in the rapid proton capture process (rp-process) play a critical role in determining energy 
release and final isotopic abundances in X-ray bursts 0j S SIS IE H- Reliable rates are therefore important for 
quantitative interpretations of observations. For example, new highly accurate data on burst profile changes over 
periods of years as observed in GS 1826-24 would provide unique constraints on X-ray burst models if the nuclear 
physics of the rp-process vifould be well enough understood. As demonstrated in a number of X-ray burst model 
calculations 0, 0, Ojj this is clearly not the case. 

One problem is the reliability of proton capture rates in the rp-process, especially in X-ray bursts where the 
reaction path runs close to the proton drip line. Far from stability, proton capture rates typically are governed by 
a few resonances and therefore statistical models are not applicable over the entire relevant temperature range of 
0.2-2 GK JUl. On the other hand, shell model calculations in the sd- [Ulill and fp-shells [H El 113 predict the 
properties of individual states but the typical uncertainty in level energies is about 100 keV even for the calculations 
of level shifts from mirror states. This translates into uncertainties of 3 or more orders of magnitude in the reaction 
rates [ll Il8|. Currently, we expect that the majority of the proton capture rates in the rp-proces suffer from such 
uncertainties. 

Clement et al. recently developed a new experimental method using radioactive beams at the National Super- 
conducting Cyclotron Laboratory at Michigan State University to accurately determine excitation energies of nuclei 
in the rp-process path. They presented new results for the excitation energies of states in ^'^Ar that determine the 
'^^Cl(p,7)'^^Ar reaction rate. Based on the new experimental data we present here a reevaluation of the •^2Cl(p,7)^'^Ar 
rate, which is part of the reaction flow through the ^°S - ■^^Ar region, an important bottleneck in the rp-process in 
X-ray bursts possibly related to the observed double peaked structure of some X-ray burst profiles . 

The '^2Cl(p,7)'^^Ar reaction rate recommended in reaction rate compilations ^3 ^ind used in X-ray burst models 
so far was a ground state capture rate based on shell model calculations with some experimental information from 
the 33 Ar mirror nucleus ^^p jl4j . We use the new experimental data on states in ^^Ar together with shell model 
calculations and a direct capture model to derive a new dramatically improved rate for the ^^C\{p,j)^^At reaction. 
This differs from the initial discussion in Clement et al. as we also take into account the contribution from proton 
capture on the first excited state in ^^Cl at 89.9 keV. This state is thermally populated in an astrophysical plasma and 
as we will show here actually dominates the proton capture rate on ^^Cl for the important temperature range between 
0.2 and 1.5 GK. We also implement slight improvements in the penetrability calculations, evaluate the remaining shell 
model uncertainties and present the reaction rate in a form usable by astrophysical reaction networks. 
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TABLE L Properties of resonant states. Listed are spin and parity J'^ , excitation energy Ex, center of mass resonance energy 
Er, proton single particle widths Fsp for angular momenta I, spectroscopic factors C^S, proton-decay width Fp, 7-decay width 
F-y and the resonance strength C1J7. The upper part is for ground state capture, the lower part for capture on the first excited 
state in ^^Cl. 



J" 


En-iMeV) 


Er{MeV) 


r 

A sp 

1=0 


1=2 


1=0 


s 

1=2 








5/2+ 


3.364 


0.021 ± 0.009 




8.00e-42 




2.90e-02 


2.32e-43 


1.77e-02 


2.32e-43 


7/2+ 


3.456 


0.113 ± 0.009 




1.90e-13 




3.00e-03 


5.70e-16 


1.94e-03 


7.60e-16 


5/2+ 


3.819 


0.476 ± 0.008 




2.10e-02 




4.30e-02 


9.03e-04 


1.54e-02 


6.88e-04 


1/2+ 


4.190 


0.847 ± 0.100 


6.30e+02 


l.Ole+01 


6.70e-02 


3.40e-02 


4.26e+01 


1.54e-01 


5.11e-02 


3/2+ 


4.730 


1.387 ± 0.100 


2.50e+04 


6.30e+02 


3.00e-03 


3.90e-02 


9.96e+01 


8.48e-02 


4.33e-03 


7/2+ 


3.456 


0.023 ± 0.009 




1.30e-39 




4.64e-03 


6.03e-42 


1.94e-03 


4.83e-42 


5/2+ 


3.819 


0.386 ± 0.008 


1.35e-01 


1.50e-03 


2.40e-02 


4.39e-01 


3.90e-03 


1.54e-02 


1.78e-03 


1/2+ 


4.190 


0.757 ± 0.100 




3.50e+00 




2.27e-03 


7.94e-03 


1.54e-01 


5.73e-06 


3/2+ 


4.730 


1.297 ± 0.100 


1.60e+04 


3.80e+02 


7.49e-02 


3.77e-03 


1.20e+03 


8.48e-02 


3.13e-02 



II. RESONANT CAPTURE 



The resonant reaction rate for capture on a nucleus in an initial state i, Na < (Jv >rcsi can for isolated narrow 
resonances be calculated as a sum over all relevant compound nucleus states j above the proton threshold ^] : 

Na < (TV >rcsi= 1.540 X 10"(/^T9)"^/2 ^ w7y e--^'^/('=^) cm^s'^mole"^ (1) 

j 

where the resonance energy in the center of mass system Eij — Ej — Q — Et is calculated from the excitation energies 
of the initial Ei and compound nucleus Ej state, k is the Boltzmann constant. We adopt an experimental ground 
state mass difference Q (reaction Q- value) of 3.343±0.007 MeV [23|- Tg is the temperature in GK and fi is the reduced 
mass of the entrance channel in amu. The resonance strengths uj^ij are in MeV and can be calculated for proton 
capture as 

9 7 . _|_ 1 r . T • 

^ ^^j + i ip^jlTJ (2) 
2(2 Ji + 1) Ftotalj 

where Ji is the target spin, and Jj, Fpy, T^j, Ftotaij are spin, proton decay width, 7-decay width, and total width 
of the compound nucleus state j. The total width is given by Ftotaij = ^-yj + J2i^pij- The proton decay widths 
depend exponentially on the resonance energy and can be calculated from the proton spectroscopic factor C^Sij and 
the single particle proton widt h Ts -nn as Tpij = C'^SijTgpij- Here we calculated spectroscopic factors with the USD 
shell model as in Herndl et al. UJ. The single particle proton widths for most states were calculated from an exact 
evaluation of the proton scattering cross section from a Woods-Saxon potential well. This method is more precise 
than the penetrability model used by Herndl et al. and Clement et al. 0, but agrees with previous work within 
20%. The 7 widths F-^ were taken from Herndl et al. This is justified as new experimental level energies are 

only available for the lower lying resonances where lo^ is largely independent of F^ as F-y >> Fp, and because for 
those states the changes in the 7-widths are less than 20%. The resulting properties of the resonance states in '^'^Ar 
are listed in Tab. 



III. DIRECT CAPTURE 



The contribution to the '^^Cl(p,7)"^'^Ar rate from direct capture into bound states has been calculated with a potential 
model using a Woods-Saxon nuclear potential (central-|-spin-orbit) and a Coulomb potential of a uniform charge 
distribution. The nuclear potential parameters were determined by matching the bound state energies. Spectroscopic 
factors were calculated with the USD shell model as in Herndl et al. [T^ . 

The direct proton capture capture rate Na < crv >dci on the target nucleus in state i is usually parametrized in 
terms of the astrophysical S-factor S{Eo) in MeV barn at the Gamov window energy Eq p^ : 

Na < av >dci= 7.83 x 10^ ( S,{Eo)e-^-^'>^^"^/^'^"'" cm^ s" Vole"! (3) 
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TABLE II: Spectroscopic factors C^S and astrophysical S-factors S{Eo) for direct capture into bound states in '^''Ar. Listed 
are results for capture on the ^^Cl ground state as well as on the first excited state in ^^Cl (denoted with an asterix). are 
spin and parity of the '^''Ar final state, n is the node number, lo the single particle orbital momentum, and jo the total single 
particle angular momentum. The table only includes the significant contributions to the total S-factor. 





J" 


(n/o)jo 


C^S 


5'(£o)(MeVbarn) 




5'(£o)*(MeVbarn) 





1/2+ 


2Sl/2 


0.08 


6.56e-3 










ld3/2 


0.67 


4.47e-3 


1.13 


4.45e-3 


1.359 


3/2+ 


2Sl/2 






0.006 


6.61e-4 






ld3/2 


0.19 


1.78e-3 


0.12 


6.80e-4 


1.798 


5/2+ 


2Si/2 






0.002 


3.21e-4 






ld3/2 


0.15 


1.84e-3 


0.62 


4.72e-3 






ld5/2 






0.021 


1.65e-4 


2.439 


3/2+ 


2Sl/2 


0.031 


4.32e-3 


0.024 


1.99e-3 






ld3/2 


0.17 


1.09e-3 


0.13 


5.01e-4 


3.154 


3/2+ 


2Si/2 


0.068 


1.02e-2 










ld3/2 


0.52 


2.21e-3 


0.17 


4.35e-4 


total 








3.25e-2 




1.39e-2 



with Z being the charge number of the target nucleus. Si{Eo) is the sum of the individual S-factors of the transitions 
from the initial state i into all bound states in the final nucleus. Table ITU lists the individual S-factors found for 
transitions into the 5 bound states of '^'^Ar from both, the ■^^Cl ground state and the first excited state as well as the 
total S-factor. For ground state capture, our total S-factor agrees within 30% with the one derived by Herndl et al. 
[T^ using the same spectroscopic factors. 



IV. NEW REACTION RATE 

The relative population of the target nucleus states in thermodynamic equilibrium is simply given by the Saha 
equation. The total reaction rate is then the sum of the capture rate on all thermally excited states in the target 
nucleus weighted with their individual population factors: 

(2 J -I- lle"'^"/''^ 

Na < (TV >^ 22^Na < av >resz +Na < <TV >dc^) ^^ , \ ' . _^ /fc^ • (4) 

To obtain the total resonant reaction rate on a thermally excited target one can combinine Eqs. and 2| and finds 

m 

NA<av>^Y.N^<^->o,7^M+ E^r^l (5) 

where j sums over all resonances in the compound nucleus, i sums over the thermally populated states in the target 
nucleus as long as the resonance energy Eij > with i — being the ground state. Na < crv >oj is the reaction rate 
contribution from ground state capture via the resonance j. G{T) is the temperature T dependent partition function 
of the target nucleus: 

G(T) = (2Jo + ^(2J, + l)eM-Er/kT) (6) 

i 

This is similar to Eq. 10 in Fowler, Caughlan and Zimmerman p3 |. Eq.Oshows that for a given resonance the relative 
contribution from capture on each excited target state depends only weakly on temperature and the actual population 
of the excited state through the partition function G(T). The contribution of excited states is mostly determined by 
the ratio of the proton widths to the excited state and to the ground state. Of course the temperature determines 
through < (TV >oj the relative importance of the various resonances. 

In this work we consider only capture on the ground state and the first excited state in ^^Cl. The ^^Cl ground 
state is experimentally known to have spin and parity 1+. The spin for the experimentally known first excited state 
at 89.9 keV has not been determined unambiguously but we assign a spin of 2+ based on the level structure of the 
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•^^P mirror and our shell model calculations. Excited states in the target will only play a role when for the relevant 
temperatures the thermal excitation timescale is smaller than the proton capture timescale. We can estimate the 
thermal excitation timescale for a 7 induced transition from the ground state to the first excited 2+ state in ^^Cl by 
using the formalism of Ward and Fowler ^2^. We assume a level lifetime of 280 ps from the mirror level in '^^P, which 
is an upper limit as the mirror state is 12 keV lower in excitation energy. We find that even for the lowest relevant 
temperatures of 0.2 GK the excitation timescale is of the order of 30 ns. Assuming typical X-ray burst conditions 
with a density of 10^ g/cm^ and a proton mass fraction of 0.7 the timescale for thermal excitation is always less than 
0.2% of the proton capture timescale for all temperatures between 0.2 GK and 2 GK. Therefore, the first excited 2+ 
state in '^^Cl is always in thermal equilibrium with the ground state and Eq.jS] applies. 

What is the role of higher lying excited states in ^^Cl? The second excited state in ^^Cl is located at an energy of 
466.1 keV. At such a high excitation energy the proton decay width Fp in Eq. [Sjis already drastically reduced for the 
low lying resonances that have a significant ground state capture rate < av >oj- Therefore, the capture rate on the 
second excited state in ^^Cl is negliglible as either Fp2j or < av >oj in Eq. (Sjare small for all resonances. Similar 
arguments apply to higher lying states. 

The contributions from individual resonances and direct capture to the total '^^Cl(p,7)'^'^Ar reaction rate are shown 
in Fig.^ The temperature dependent relative population of the ground and first excited state in ^^Cl has been taken 
into account. Direct capture on the ground or the excited state is negligible over the entire relevant temperature 
range of 0.2-2 GK. 




temperature (GK) 



FIG. 1: The contributions of various individual resonances and through direct capture (DC) to the ^^Cl(p,7)^^ Ar reaction 
rate as functions of temperature. In the legend, resonances are labeled with their excitation energy in "^"^Ar. The upper panel 
shows contributions from ground state capture, the lower panel contributions from capture on the first excited state in ^'^Cl in 
each case weighted with the relative population of the respective target state. In addition, both panels show the same total 
^^Cl(p,7)'''^Ar reaction rate for comparison. Contributions from ground state capture via the 3.36 MeV resonance and from 
capture on the first excited state via the 3.46 MeV resonance are too small to appear in this graph. 

As Fig. n shows, the capture rate contributions from ground and excited target states vary greatly as Fp depends 
strongly on proton energy, spin, and nuclear structure. For '^■^Cl(p,7)'^^Ar the most dramatic change occurs for the 
5/2+ 3.819 MeV resonance. For this resonance the spin of the excited state (2+) is one unit larger than for the ground 
state (1+) allowing for s-wave protons to populate the 5/2+ resonance in addition to d-wave protons. As a result the 
proton width for capture on the excited state is a factor of 4.3 larger than for ground state capture despite the slightly 
lower proton energy. As a consequence (see Eq. ISJ the resonant capture on the first excited state in ■^^Cl via the 3.819 
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TABLE III: Spectroscopic factors C^S for neutron removal from '^^Ar to various 5/2^ states in '^'^Ar calculated with the 
shell model effective interactions USD and USD05 and compared to experimental data 28] from the '^'*S to '^''P mirror reaction 
for d5/2 protons. The excitation energies Ex of the 5/2^ states are the experimentally known ones in the ^^P mirror. 



Ex (MeV) C^S C^S C'^S 

Exp Exp USD05 USD 

1.85 1.09 0.99 1.09 

3.49 0.31 0.27 0.70 

4.05 1.28 1.62 1.30 

5.05 1.66 1.51 1.33 



MeV 5/2+ state in '^^Ar becomes the dominant contribution to the stellar '^^Cl(p,7)^'^Ar reaction rate for all relevant 
temperatures between 0.15 and 2 GK (see Fig. This is despite the fact that the population of the first excited 
state in '^^Cl for example at 0.3 GK is only 5%. As Eq. |S1 demonstrates, the relative importance of contributions 
from capture on various states in the target is largely independent of their population (with the exception of their 
contribution to the partition function), what matters most is the proton width for the transition into the resonant 
state. 

The remaining uncertainty from resonance energies ranges from 20 - 50% for the most important temperatures 
between 0.3 and 0.8 GK and increases to about a factor of 2 at 1.5 GK. It originates mainly from the small experimental 
uncertainties in the ^'^Ar excitation energies and the reaction Q-value. For the experimentally unknown higher 
lying levels we assume an uncertainty in the excitation energy of 100 keV, which is a typical deviation for shell 
model predictions of Coulomb shifts in this region. These levels only begin to contribute to the total reaction rate 
beyond around 1 GK and lead to the somewhat increased error at highest temperatures. At this level of error, other 
uncertainties in the shell model calculations now become relevant as well. The Woods-Saxon well used for the proton 
scattering widths is adjusted so that the total density for protons filled up to the Fermi energy reproduces the rms 
charge radius of nearby stable nuclei. We estimate an uncertainty of 5% in the Woods-Saxon radius paramter which 
leads to a 15% uncertainty in the proton decay widths. 

In principle one could estimate the uncertainty of the shell model spectroscopic factors by comparing calculated 
spectroscopic factors with experimental data. No experimental data are available for the here relevant proton spec- 
troscopic factors in '^'^Ar or the corresponding neutron spectroscopic factors in the ■^'^P mirror. In addition, small 
spectroscopic factors are in general difficult to reliably extract from direct reaction cross sections due to the possibil- 
ity of multi-step routes. However, we can test our shell model descriptions for the relevant states in ^'^Ar by comparing 
their relatively large neutron spectroscopic factors for '^'^Ar to ^'^Ar neutron removal with experimental data for the 
mirror proton spectroscopic factors for ^''S to ^'^P proton removal. We limit the discussion to the 5/2"'" states, as the 
■^^Ar(p,7)'^^Ar is dominated by the resonant contribution from the third 5/2+ state in ^'^Ar. The calculations and 
experimental data for the 5/2+ states are summarized in Tab. IIIll The calculations are carried out with the original 
USD interaction as well as close to final version of a new sd-shell interaction USD05 |23| obtained from a least squares 
fit to about 600 levels in the sd-shell nuclei. The spectroscopic factors are all large numbers and in good agreement 
with the experimental data from the mirror reaction. When calculated levels can be matched to experimental levels 
with the present degree of accuracy (about 100 keV) the comparison in Tab. IIIII shows that the largest calculated 
spectroscopic factors are accurate to about 20% and the moderately large spectroscopic factors (as for the second 
5/2+ state) are accurate to about a factor of two. 

In addition, we can look at variations in spectroscopic factors calculated with different shell model effective inter- 
actions. The spectroscopic factors for the important ^^Cl (2+) to ^'^P (5/2+) I = transitions are shown in Tab. II VI 
for the old USD interaction as well as the new USD05 interaction. All of the spectroscopic factors are small since 
the Si/2 orbit is mostly filled in ^^Cl. The 1 = spectroscopic factor summed over all final 5/2+ states is only 0.131 
for USD05. The comparison of the two interactions in Table 2 shows that the largest of the spectroscopic factors 
(including the one to the third 5/2+ state of interest) agree to about 20 %. Spectroscopic factors calculated with 
other effective interactions such as SDPOTA or CWH [33|, which also do a reasonable job of reproducing the 
energy levels for nuclei in the upper sd-shell, show a similar behavior. In general a small spectrosopic factor in the 
presence of other states with much larger spectroscopic factors may be very uncertain (as observed for many of the 
states in Tab. II V|) . However, when the spectroscopic factor is large compared to those for nearby states it can be 
relatively accurate. Thus the estimated theoretical error for small spectroscopic factors must be treated on a case by 
case basis. 

To summarize, the shell model can describe the neutron single particle strength of the here relevant states within 
20% to a factor of 2. In addition, the important I = proton spectroscopic factor for the third 5/2+ state is relatively 
large compared to other states and the total / = proton single particle strength, and its calculation should therefore 
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TABLE IV: I — spectroscopic factors for the 5/2'^ states in ''"^Ar for proton capture on the 2'^ state in "'^Cl. Spectroscopic 
factors C^S and excitation energies Ex calculated with the USD and the USD05 interaction are compared. 



Ex (MeV) 


Ex (MeV) 


C'^S 


C'S 


USD05 


USD 


USD05 


USD 


1.99 


2.00 


0.0013 


0.0006 


3.42 


3.83 


0.0022 


0.012 


4.22 


4.19 


0.025 


0.023 


4.82 


5.00 


0.049 


0.040 


6.39 


6.64 


0.00002 


0.0024 


6.54 


6.93 


0.017 


0.0024 



TABLE V: Reaction rate Na < av > as a, function of temperature T. Given is the recommended rate from this work (rec) as 
well as a lower and upper rate that reflect the estimated error bar. 



T (GK) NA<av> (cm^s/mole) 





rec 


lower 


upper 


0.1 


7.85e-15 


1.73e-15 


2.81e-14 


0.2 


6.62e-09 


3.20e-09 


1.36e-08 


0.3 


3.44e-05 


1.94e-05 


6.04e-05 


0.4 


2.09e-03 


1.27e-03 


3.41e-03 


0.5 


2.21e-02 


1.33e-02 


3.67e-02 


0.6 


9.97e-02 


5.80e-02 


1.77e-01 


0.7 


2.82e-01 


1.58e-01 


5.37e-01 


0.8 


6.06e-01 


3.27e-01 


1.24e+00 


0.9 


1.09e+00 


5.67e-01 


2.39e+00 


1.0 


1.76e+00 


8.75e-01 


4.06e+00 


1.5 


7.94e+00 


3.40e+00 


2.03e+01 


2.0 


1.78e+01 


7.48e+00 


4.39C+01 



be relatively reliable. We therefore estimate a factor of 2 uncertainty in the shell model spectroscopic factor in this 
specific case. To take into account other uncertainties such as the proton single particle widths discussed above, 
we adopt a total uncertainty of a factor of 2.5 for the ■^^Cl(p,7)'^'^Ar reaction rate on top of the uncertainty from 
the resonance energies, which is accurately calculated as a function of temperature. The new reaction rate with its 
uncertainties is tabulated in Tab. and displayed in Fig. [21 The total estimated uncertainty of the rate is about a 
factor of 3 at temperatures around 0.3-0.7 GK to about a factor of 6 at temperatures around 1.5 GK. 




temperature (GK) 



FIG. 2: Total Cl(p,7) Ar reaction rate from this work (solid area) and previous shell model rate without experimental data 
(hatched area). The width of the area reflects the estimated uncertainty. 
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TABLE VI: Fit coefBcents (see Eq. |7| for the recommended rate from this work (rec) as weU as for the lower and upper limits 
that reflect the estimated error bar. 



rec 


ao 

0.150730E+01 
0.542869E+02 


11 

-0.539842E+01 
-0.782163E+01 


0.2 

-0.153745E+02 
0.367513E+03 


03 

0.187918E+02 
-0.425108E+03 


14 

0.126139E+01 
0.502875E+01 


14 

-0.208373E+00 
0.823457E+00 


17 

-0.134957E+02 
0.262097E+03 


lower 


0.551123E+03 
-0.517996E+02 


-0.158099E+02 
-0.156901E+01 


0.886227E+03 
-0.173733E+03 


-0.150891E+04 
0.237842E+03 


0.878013E+02 
-0.114259E+02 


-0.495734E+01 
0.538535E+00 


0.714263E+03 
-0.124462E+03 


upper 


0.299493E+03 
-0.269317E+03 


-0.143477E+02 
-0.794584E+00 


0.705984E+03 
-0.279531E+03 


-0.104556E+04 
0.597651E+03 


0.530983E+02 
-0.505961E+02 


-0.269508E+01 
0.399204E+01 


0.533394E+03 
-0.244496E+03 



Often, reaction rates are implemented into reaction network codes using a fit of the form: 

Na < crv >= ^ exp(ao.t + auTg^ + a2iTg^^^ + a^iTg^^ + a^Tg + a^iT^^^ + aeJnTg) (7) 

i 

which is for example used in the reaclib reaction rate library. In Tab. IVII we give fit parameters a„i for our new 
recommended •^^Cl(p,7)'^^Ar rate as well as for the upper and lower uncertainties. The fits are accurate to 5% and 
are vahd between 0.1 and 10 GK. Beyond 2 GK we use the Hauser-Feshbach predictions from NON-SMOKER 
using the procedure described by the NACRE collaboration [sj- The low temperature behavior is reasonable, so the 
fits can be used in model simulations that encounter lower temperatures. 

V. DISCUSSION 

Our ground state capture rate for ■^^Cl(p,7)'^^Ar is similar to the rate given by Clement et al. |l|. Small differences 
arise from the slightly improved calculation of the proton single particle widths, the comprehensive reevaluation of 
the shell model uncertainties, and the modifications in the proton widths by taking into account transitions into 
excited target states. The major change however comes from the consideration of the capture on the first excited 
state in '^^Cl. In the critical temperature range around 0.2-0.4 GK where for typical X-ray burst conditions (densities 
of 10^^^ g/cm'^ and hydrogen mass fractions of 0.1-0.7) the reaction rate becomes comparable to the burst timescales 
(10-100 s) our rate is a factor of more than 4 larger than the ground state capture rate shown in Clement et al. 0]. 
In the same temperature range our rate is up to a factor of 70 larger than the recommended shell model based rate 
in Herndl et al. |l| (see Fig. Ej). 




temperature (GK) 

FIG. 3: Ratio of the ''^Cl(p,7)'^''Ar reaction rate of this work to the rate calculated by Herndl et al. [T4|. 

The impact of thermally excited states in the target nucleus on a reaction rate is often expressed in terms of a 
"Stellar Enhancement Factor", or SEF, which is defined as the ratio of the actual capture rate to the ground state 
capture rate. Fig. ^ shows the SEF determined in this work. The enhanced proton capture rate from the first 
excited state in '^^Cl leads to a dramatic enhancement of the reaction rate of up to a factor of 5 at the important 
temperatures around 0.2-0.4 GK. Recent reaction rate compilations typically use SEFs that are calculated with the 
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FIG. 4: Stellar enhancement factor derived in this work (solid line) and from the statistical model code NON-SMOKER 
(dashed line). 

statistical Hauser-Feshbach mod el fl9ll3l| . For comparison we also show in Fig.^jthe SEF obtained with the Hauser- 
Feshbach code NON-SMOKER Clearly in this case a statistical approach to calculating the SEF does not give 
the correct result, even in the temperature range above 0.7 GK where strictly speaking the Hauser-Feshbach model 
is expected to be applicable [13, |33 ■ This is no surprise as here the SEF is entirely dominated by the properties of a 
single resonance. However, this example illustrates that large uncertainties can be introduced when using statistical 
model SEFs far from stability when reaction rates are dominated by a few resonances only. A better general approach 
would be to calculate SEFs within the shell model. 

Clearly accurate masses and excitation energies are the single most important step for improving the accuracy of 
theoretical reaction rates. For comparison Fig. |21 also shows the uncertainty band of the reaction rate for the case that 
none of the excitation energies would be known experimentally, a situation that is very common along the rp-process 
path. Clearly, without experimental excitation energies reaction rates in the rp-process far from stability can be 
uncertain by many orders of magnitude. A measurement of excitation energies and Q-values to better than 10 keV 
can reduce this uncertainty to about a factor of 3. 

To obtain in critical cases even more precise reaction rates one would need to perform indirect measurements of 
spectroscopic factors, or perform a direct measurement of the reaction rate at astrophysical energies. Note, however, 
that in this particular case a direct measurement would not be possible as the reaction rate is dominated by capture on 
the first excited state in the '^^Cl target nucleus. This underlines the importance of indirect methods in determining 
accurate stellar reaction rates. 

We thank M. Wiescher for pointing out the importance of the timescale for thermal excitation and O. Sorlin for 
pointing out the potential relevance of the first excited state in '^^Cl. This work was supported by NSF grants PHY 
0110253 (NSCL) and PHY 0216783 (Joint Institute for Nuclear Astrophysics). H.S. acknowledges support through 
the Alfred P. Sloan Foundation. B.A.B. acknowledges support from NSF grant PHY-0244453. C.A.B. thanks for 
support by the U.S. Department of Energy under grant No. DE-FG02-04ER41338. 
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